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2 YUNBO ZENG AND WEN-XIU MA 

I. Introduction. 

As is known, some integrable systems possess bi-Hamiltonian structure. We recall 
some known results. Let M be a differential manifold, TM and T*M its tangent and 
cotangent bundle, and ^o and 6i : T*M ^ TM two compatible Poisson tensors on M 
[1]. A vector field X is said to be bi-Hamiltonian (BH) with respect to 9o and 9i, if two 
smooth functions, H, F E C°°{M), exist such that 

X = OodH = OidF (1.1) 

where dF denotes the differential of F (gradient V-^ fo^^ finite system and variation 5F 
for field system). If ^o is invertible, the tensor (^ = 9i9q is a Nijenhuis tensor or hered- 
itary operator. The operator $ maps a given BH vector field into another BH vector 
field. Hence having a Nijenhuis tensor, one can construct a hierarchy of Hamiltonian 
symmetries, and a related hierarchy of integrals of motion for the underlying system. 
The BH structure (1.1) ensures that the resulting integrals of motion are pairwise in 
involution with respect to both Poisson brackets. Thus the BH structure of a given 
system is important for its integrability. 

Unfortunately, for a majority of the BH finite-dimensional systems, none of the ^o 
and 6i is invertible. In fact, all the known BH finite-dimensional systems arising from 
the constrained fiows or stationary fiows of soliton equations usually exist in an extended 
phase space and both ^o and 6i are degenerated (see, for example, [2-8]). In their natural 
phase space these systems may satisfy a weaker condition than the BH one. The notion 
of a quasi-bi-Hamiltonian (QBH) system was introduced [9,10]. According to [10], for 
dimM = 2n, a vector field, X, is said to be a QBH vector field with respect to Poisson 
tensors, ^o and ^i, if there exist three smooth functions H, F, p, such that 

X = ^oV^=-^iV^ (1-2) 

P 

where two Poisson tensors ^o and 6i are compatible and nondegenerated (invertible). 
The function p is called an integrating factor. On a 2n-dimensional symplectic manifold 
M, let (g = {qi, ■■■,qn),P = {pi, ■■■^Pn)) be a set of canonical coordinates and Oq the 

canonical Poisson matrix ^o = I t n ] (^ denoting the n x n identity matrix). As 

6*0 and 9i are compatible and invertible, the Nijenhuis tensor $ = 9i9q is maximal, 
i.e. it has n distinct eigenvalues /x = (/xi, ...,fXn)- As is known [11], in a neighborhood 
of a regular point, where the eigenvalues /x are distinct, one can construct a canonical 
transformation {q,p) i— > (/x, t')((/i, i^) referred to as the Nijenhuis coordinates) such that 
9i and $ take the Darboux form 

^i=f ^^ qM , ^=(o^ aJ' ^i = diag{pi,...,pn). (1.3) 
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A QBH vector field is said to be Pfaffian [10] if, in the Nijenhuis coordinates, an inte- 
grating factor p in equation (1.2) is the product of the eigenvalues of $, i.e. 



n 



Hi- (1.4) 



In the Pfaffian case, the general solutions, H and F, of equation (1.2) are obtained 
and the Hamilton- Jacobi equation for H is shown to be separable by verifying the Levi- 
Civita conditions [12]. Some relationship between BH and QBH structure is discussed 
in [13]. Several QBH systems are presented [9,10,12-14]. It is in general quite difficult 
to directly construct a BH or QBH structure for a given integrable Hamiltonian vector 
field. In recent years much work has been devoted to the constrained fiows of soliton 
equations (see, for example, [2-8,15-24]). One of the aims of this paper is to show how 
to construct an infinite number of families of QBH systems from the constrained fiows 
of soliton equations. We have presented some families of the constrained fiows in order 
to study the dynamical r-matrices in [24]. We now describe the explicit QBH structures 
for these families of the constrained fiows. 

The Lax representation for the constrained fiows of soliton equations can always be 
deduced from the adjoint representation of the Lax pair for soliton equations [16,17]. 
There is an effective way for the separation of variables for some finite-dimensional in- 
tegrable Hamiltonian systems with some kind of Lax matrices [25,26]. The separated 
variables for some constrained fiows can be introduced and the Jacobi inversion prob- 
lems for the constrained fiows can be established by means of the Lax representation 
[27,28]. We are interested in the relationship between the two methods for the separa- 
bility mentioned above. Another main aim of this paper is to prove that the Nijenhuis 
coordinates for the underlying families of QBH systems are usually the same as the 
separated variables introduced by the Lax matrices. 

The paper is organized as follows. 

In section 2 we present a new QBH system. We directly construct the second compat- 
ible Poisson tensor by using a map relating this system to its modified version, and prove 
that the Nijenhuis coordinates for this system is equivalent to the separated variables 
defined by Lax matrix. We make some comparison of the two methods for separability. 
In section 3 and section 4, by using the constrained flows associated with the polynomial 
second order spectral problems and the higher-order symmetry constraints, we propose 
a way to construct an infinite number of families of QBH systems. The explicit QBH 
structures of the first two families of constrained fiows are given. The equivalence of the 
Nijenhuis coordinates and the separated variables is proved. In section 5 we point out 
that the two compatible Poisson tensors 6*0, Oi and the integrating factor p given by the 
QBH structure (2.28) and (2.29a) are just that for the third family of OBH systems. 
Also some conclusions and a conjucture are given. 

II. New QBH system. 
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In this section we present a new QBH system. By using a map relating this system to 
its modified version, the second compatible Poisson tensor is obtained from the image 
of the Poisson tensor for the modified version under the map. We use this system to 
illustrate how to prove the equivalence of the Nijenhuis coordinates and the separated 
variables introduced by the Lax matrix. 

A. New finite-dimensional integrable Hamiltonian system. 

For Jaulent-Miodek (JM) spectral problem [29] 

^. =£/(„, AW,, ^(«. A) =(,._„»,_„„ I), ^-(^:), «=(:;). (2.1) 

its adjoint representation is defined by [30] 

V^ = [U,V] = UV-VU, (2.2) 

where V is taken as 

tti hi 



i=0 



V = J:V.X-'. V.= r; _"• . (2.3) 



Then equation (2.2) and (2.3) yields 



ao = ai= a2 = bo = bi = 0, 62 = 1, ^3=2"!' 
1 1 

«3 = -tMi,x, Co = 1, Ci = --Wi,..., 



and in general 

^bk+ij \ h J 



where 



(2.4a) 



ak = ~-^bk,x, Ck = ak,x - uobk - uibk+i + bk+2, k = l,2,---, (2.4b) 



\ 1 / ' 9x' 

The Jaulent-Miodek hierarchy associated with (2.1) can be written as an infinite- 
dimensional Hamiltonian system 
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where the Hamiltonian H^ and the Hamiltonian operator J are given by 

•^ = (2D -uJ^2u,d) ' ^- = n^^^-+' " ''^^-+'^- 

Under zero boundary condition we have 

The constrained flow of (2.5) consists of the equations obtained from the spectral 
problem (2.1) for A^ distinct \j and the restriction of the variational derivatives for 
conserved quantities Hi (for any fixed /) and Xj [15-17]: 

^i,x = ^2, ^2,x = A2*i-'UiA^i-tio^i, (2.7a) 

5Hi 1 A 5\. fbi+2\ 1 /< A*i, *i > 



E^- r -^r;:.";rho. (2.7b) 



5u 2 ^ 5u \bi+ij 2V<^i,^i> 

which has been recognized as a symmetry constraint [18-20]. Hereafter we denote the 
inner product in R^ by < ., . > and 

^i = ii^iir-- ,i^iNV, i = 1,2,A = diag{Xi,--- ,Xn)- 
For / = 4, we have 

^^ ^ 128^^ ^ 16^^^° ~ 32^^"^^^ ^ 8^°^^ ~ g^^ix^ox- (2.8) 

By introducing the Jacobi-Ostrogradsky coordinates 

gi =ui, q2 =Uo, 

Pi = T = -zr^uiuix - -uox, P2 = T = --Mix, (2.9) 

ouix Id 8 ouqx 8 

the equations (2.7) for / — 4 are transformed into a finite- dimensional Hamiltonian 
system (FDHS) 

ajp fjTp arp 

^1- = «^ = ^2, qix = ^ = -8p2, q2x = ^ = -Spi + 20qip2, (2.10a) 

0^2 Opi Op2 

arp 

^2x = -^ = A2^i-giA^i-Q2^i, (2.10b) 

aFi 35 4 15 ^ ^3^1. , , 

Pix = —g^ = Y^qf + Y^qh2 - lOpl + -ql --< K^u'^^ >, (2.10c) 
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dFi 53,3 1 



P2x = --Q-^ = Y^qt + -^qiq2 - - < *i, *i >, (2.10d) 

or equivalently 

where 

p_ (xJ.T n n iXl'^ rn r. \T Q - ( ^ ^{N+2)x{N+2) 

^ — [^i,qi,q2,^2iPi^P'2) , ^0 — r n 

\~-^{N+2)x{N+2) U 

111 1 

-8pip2 + lOqipj - —qfq2 - -qiql - —ql (2.11) 

The Lax representation for FDHS (2.10) can be deduced from the adjoint representation 
(2.2) by using the method in [16,17] which is sketched as foUows. Due to (2.4a), (2.6) 
and (2.7b), we may define 

^m = ^<A"^-'^i,^i>, m = 5,6,..., 

which together with (2.4b) and (2.10) yields 

_ 1 _ 1 

am = -2 < ^"""^^1,^2 >, c^ = -- < A"'"^*2,^2 >, m-5,6,.... 

Set 

Then the construction of a^, bm,Cm ensures that under (2.10) 

also satisfies (2.2). Notice that 



set 



FAMILIES OF QUASI-BI-HAMILTONIAN SYSTEMS AND SEPARABILITY 7 

we have 

^(A) = 2p2A + 2pi - 2gip2 - ^ f] y^, (2.12b) 



C(A) = A^ - igiA^ - (^52 + ^Q?)A2 + (^g? + ^^1^2 - ^ < ^1, ^1 >)A 

Af / 2 

+ i^2 - ^'^t -4pi - ^ < A*i,*i > +lqi < *i,^i > -IJ2yZY.- (2.12d) 

Since V under (2.10) satisfies (2.2), then Q under (2.10) satisfies (2.2), too, namely 

Qx = [t/,Q], (2.13) 

which presents the Lax representation for (2.10). This can also be verified by a direct 
calculation. The equation (2.13) implies that ^TrQ'^iX) = A'^{X) + B{X)C{X) is the 
generating function of the integrals of motion for (2.10). We have 

A\X) + B{X)C{X) = A^ - FiA + F2 + V -, (2.14) 

111 

^ Zj ^ 

+{pf + -^1iQ2 - - < ^1, ^1 >) < ^1, ^1 > --qi < ^2, ^2 > 

Q -1 Fi 1 "1 "1 

+(g^i + ^q2){-^pl - ^qf + -^ql - - < A^i, ^1 > +-Q1 < *i, ^1 >) 

1 1 

-(7^2 + —ql) < A^i, ^1 > -2p2 < ^1, ^2 > +4(pi - glP2)^ (2.15) 

4 lb 

F« = (-2P2A, + 2gip2 - 2pi)V^H^2z - -(A? + -QiA, + -ql + -§2)^1, 
+ 2 [^^ - 2^1^? " ^8^' ^ 2^'^^^' ^ ^4^' + 2^'^' - - < ^1, ^1 >)A, + -ql - Apj 

(2.16) 
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where F^^\ i = 1, ..., A^, Fi, F2 are A^ + 2 independent integrals of motion for (2.10). By 
means of the r-matrix, it can be shown that the equation (2.10) is a finite- dimensional 
integrable Hamiltonian system (FDIHS). 

In order to find the QBH structure for (2.10), we need to use the modified system of 
(2.10). Let us consider the modified Jaulent-Miodek (MJM) spectral problem [31] 

,.^Ui.,m Ui.A)^{,\ 4). ^^(2). .^C;). (2.17) 
The equations (2.2) and (2.3) yield 

ao = 0, 60 = 1, bi = -vi, ai = vo, cq = 1, ci = --vi,... 

^ ^f -2t>o + D 

\\D + ^D-^voD ^vi + ^D-^viD 

The MJM hierarchy associated with (2.17) can also be written as a infinite-dimensional 
Hamiltonian system 

where the Hamiltonian H^ and the Hamiltonian operator J are given by 

-^=Mq -2D ) ' ^" ^ — [an,x - vibn + 2bn+i]. 

Also we have 

In the similar way as for (2.7), the constrained fiow of (2.19) is defined by 

$i,x = t'o$i + A$2, $2,x = (A- vi)$i -t^o$2, (2.21a) 

5Hi 1 /2 < $i,$o >\ 

^ + - / k =0, (2.21b 

where $i = (0ii, ..., 0iAr)'^,i = 1,2. 
For / = 3, 

TT / ^ 2 2,4,^4 '^ 2 '^2 2\ 
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By introducing the Jacobi-Ostrogradsky coordinates 

qi=vi, q2 = vo, (2.22a) 

SHs 1 ^ 6Hs 1 3 2 ...... 

Pi = -T = -g^ix, P2 = -T = 77^0x-gi^i, (2.22b) 

ovix 8 dvox 2 8 

the equations (2.21) for / — 3 are transformed into a FDHS 



*i.- = «^ = §2*1+ A$2, g'lx = ^ = -8pi, 52x = Tp^ = 2p2 + Tg?, (2.23a) 

a$2 c/pi t^P2 4 



or 



where 



$2,x = -^=A$i-gi$i-§2$2, (2.23b) 

Pix = --^ = -^qiP2 - piqi - -g? - - < $1, $1 >, (2.23c) 

P2x = -^ = ql- Tfe- < $i,$2 >, (2.23d) 

dq'2 4 

Px=^0V^l, 

P- ($f,gi,g2,$^,Pi,P2)'^, 

+§2 < $1, $2 > +2 ^ A*2, $2 > "2 < A$i, $1 > +-^1 < $1, $1 > . 

B. The QBH structure for the FDIHS (2.10). 

We now estabhsh a map relating FDIHS (2.10) to (2.23), then use the map to con- 
struct the second compatible Poisson tensor for the FDIHS (2.10). 

It is known [31] that a gauge transformation between the JM and MJM spectral 
problem is as follows 

^1 =(/>!, ^2 = A(/)2 +^001, Ui=Vi, Uq = -vq^-vI, (2.24) 

which, together with (2.9) and (2.22), gives rise to the map relating (2.10) to (2.23), i.e. 
P = M{P): 

^i = $i, ^2 = A$2 + g2$i, qi = qi, 

q2 = -'2p2- -rqf-qi, pi = 5ipi + -gf + -§2P2 - 7 < $i, $2 >, P2 = Pi- (2.25) 



10 
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M'{P) 






-1^1 


q2l 





'\^l 


Pi 



(2.26) 



The map M given by (2.25) transforms all equations in (2.10) except for (2.10c) into 
the corresponding equations in (2.23) except for (2.23c). In fact, the equation (2.10c) 
with an additive constant term c = —\Fi is transformed into (2.23c) under the map 
(2.25). However, the second Poisson tensor constructed later by using the map (2.25) 
is valid for an arbitrary c, therefore we can take c = 0. The Jacobi M' of the map M 
take the form 

// OOON^ 

1 

-2§2 0-2 

$1 A 

m + \P2 -\^l qi \q2 

\ b 1 / 

According to the standard procedure [32], the image of the Poisson tensor 6q for the 
FDIHS (2.23) under the map M gives rise to the second compatible Poisson tensor for 
the FDIHS (2.10). That is 

a _ M'ft M''^ I f0iN+2)x{N+2) M 

OnxN 4^2 OatxI 

-ivl>T p2 \. (2.27b) 

OlxAf -P2 

Furthermore, by a straightforward calculation, we can show the following proposition. 
Proposition 1. The system (2.10) possesses the QBH representation 



(2.27a) 



^1=0 



A 



IxAT 



i*l 



4 
^1 



OatxI 

1 



Si 



2^f 



2^2 --g-^? -2^1 



where 



P, = eoVFi = -01 V El 
P 

3 11 

p = 5(A)U=o= gg?+ 2?2- 2 < A"'^i,^i >, 



(2.28) 
(2.29a) 



N 



El = [A\X) + S(A)C(A)]|a=o = F2-Y, K'F 



(0 



i=l 



3 1 

= (tT^Ii + 7^2)(< A-^^2, ^2 > - < A^i, ^1 >) 
lb 4 

3 1 

(tT^Ii + 7^152) < *i, *i > +2{pi - qiP2) < A-i*i, *2 > 
lb 4 

5 11 1 

+ (2p| + — gt - _g2 ^ _ < A^i, ^1 > --qi < ^1, ^1 >) < A-1^1, ^1 > 

H^ql + 2?2)(|gi - 4p2 - —ql) + 4(pi - qip2f 



+ -[< A"^^i,^2 >^ - < A"^^i,^i >< A"^^2,^2 >]• 



(2.29b) 
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C. The Nijenhuis coordinates. 

We now prove that the Nijenhuis coordinates for QBH system (2.28) are the same as 
the separated variables defined by means of the Lax matrix (2.12b). As ^o and 6i are 
compatible and invertible, the matrix 0i6q is maximal, it has A^ + 2 distinct eigenvalues 
jj, = {pi, ...,iJ.N+2)- The explicit form of the canonical transformation from P to the 
Nijenhuis coordinates (/i, u) is given in what follows. The eigenvalues ni, ...,fiN+2 are 
defined by the roots of the equation 

/(A)=| A/-^i 1=0, (2.30) 

which, since Ai depends only on (\l/i, gi, (72), gives rise to 

//, =/,(*!, gi, 92), i = l,...,iV + 2, (2.31) 

V'li = 9jilJ'), 3 = 1, ■■■,N, qi = gN+M, Q2 == gN+2{lj)- (2.32) 

Then we introduce the generating function by 

TV 

S = ^i/j2jgj{n) +pigN+M +p2gN+2{lJ'), (2.33a) 

i=i 

such that 

9S . , ,, dS dS 

0^2j OPl OP2 

dS ^^ dgj dgN+i ^ dgN+2 . ^ , . . „ /o qq n 

"^■^a;i- = g^^^a;^+^^^;^+^^^;^' ^ = i'-'^+2. (2.33c) 

The equations (2.33b) reconstruct (2.31) or (2.32), the equations (2.33c) give the expres- 
sion for Uj. The system (2.10) written in terms of (/x,!^) can be shown to be separable. 
On the other hand, the separated variables (/i, i^) for (2.10) can be constructed by 
means of the Lax matrix in the following way [27,28]. The coordinates jli, ...,Jin+2 are 
introduced by the zeros of -B(A): 

i.(A)^A^4,.A.?,?.i„.ii:^^||^. (2,34) 

with 

N+2 N+2 N N 

i?(A)= n(A-/ifc) = $^/3fcA^+'-', i^(A) = ^(A-Afc)=^a,A^-^ (2.35) 

k=l k=0 fc=l k=0 



12 
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N N 

N+2 N+2 

i=i i=i 

and the canonically conjugate coordinates i7i, ..., i'n+2 are defined by 



1 ■r;^ ll^lji^2i 

Uk = -A{fik) = -2p2/iA; - 2pi + 2giP2 + o >^ . 



/c = l,...,iV + 2. (2.36) 



The FDIHS (2.10) in terms of the coordinates {pt,i^) will be shown to be separable later. 
We have the following proposition. 

Proposition 2. The Nijenhuis coordinates (/x,!^) defined by (2.30) and (2.33) are ex- 
actly the same as the separated variables ifi^v) defined by (2.34) o'^^ (2.36). The QBH 
vector field (2.28) is Pfaffian in the Nijenhuis coordinates. 

Proof. We first show that 



/(A) = B{X)K{X) = R{X). 



(2.37) 



We denote /(A) by /Ar(A; Ai, ..., Xn) in order to prove (2.37) by induction. Obviously, 
(2.37) holds for A^ = 1. Then we have by induction 



/jv(A; Ai, ..., Aat) 



A- Ai 

A -As 





-2V^ii -2^12 
















A - Aat ji^iN 

A - gi -1 

-2^iN \q2 + f g? A + f gi 

TV 



(A - Ai)/Ar-i(A; A2, ..., A^) + Ul^ X{{X-X 



(2.38) 



k=2 



(A^ + J^iA + U + 1^2 + \f: -^mX) + l^KiX) = BiX)KiX). 



2 " 2 ^^ A - A, 

J =2 ^ 



2 A- Ai 



The equation (2.38) implies that Ai, similarly Xk,k = 2, ..., A^, is not the zero of /(A). 
Thus (2.37) indicates that /(A) and B{X) have the same zeros, i.e. /i^ = fik- 
It follows from (2.34) that 



,2 _ 2 -^(Aj) 
^^^•-V'(A,)' 



qi = 2(/3i -tti), 



(2.39a) 
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\ci, + ^ql = I < A-^v^i, vPi > +^^, (2.39b) 

/ o Z ttTV 

where the prime denotes differentiation with respect to A. The equations (2.39a) and 
(2.39b) yield 

^^ = 2f:-§^-3(/3,-a,)^ + 2^. (2.39c) 

According to (2.33a), one gets 



Notice that 



_^y^ , / 2.R(Aj) _ >^ tl)2jR{\j) _ 1 >^ ^ij^2j 

^^^k fri "^"1 K'{\,) - A. ^2i?(A,)K'(A,)(/x. - A,) " 2 ^tt ^^ - ^^• 

9 A 2i?(Aj) _A 2R{\j) 

J_^mX^^^2RiX^ l = J-«A-v,, v,,>+f ^" 



Ij 



IJ-k "—[XjK^Xj) {pk- Xj)K'{Xj) Hk ' ^/Wfc-Aj ' 

d/3N+2 _ Pn+2 d{l3i -Qi)^ _ _ d{/3i - ai) _ _ 

duk Hk onk OfJ'k 

and using B(nk) = 0, one finds from (2.33c) that z/fc, fik satisfy (2.36). Finally, it follows 
from (2.39b) that 

P - S(A)|a^o = 1^2 + y - \ < A-^vPi, vPi >= ^^ = (Zll^ Jj' (2.40) 

This completes the proof. 

D. Comparison of the two methods for separability. 

For the FDIHS with QBH structure, the separated variables, i.e. the Nijenhuis 
coordinates, can be introduced by the Nijenhuis tensor. Then the separability of the 
Hamilton-Jacobi equation for the system can be shown by varifying the Levi-Civita 
conditions. For the FDIHS with some kind of Lax representation, the separated variables 
can be found and the separability of the Hamilton Jacobi equation for the system can 
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be shown by means of the Lax representation. So far there is not an effective way to 
define separated variables for the FDIHSs with some kind of Lax matrices, such as the 
3x3 Lax matrices [22] or the Lax matrices admitting dynamical r-matrix. However, if 
the separated variables can be introduced by the Lax matrix, one can further establish 
the Jacobi inversion problem for the system by means of the Lax representation. By 
using the standard Jacobi inversion technique, the solution to the system can be found. 
We now use the Lax representation (2.12) to construct the Jacobi inversion problem 
for (2.10). Set 

W(\) ^+^ 

AHX) + B{X)C{X) = -^, W{X) = J2 P^>^\ (2-41) 

then Pi are also the integrals of motion for (2.10). By substituting (2.13) and using 
(2.14), (2.41) leads to 

Pn+6 = 1, PN+6-i = oti, i = 1, 2, 3, 4, 

^1 = --Pa^+i +05, ^2 = ^Af -oi^Af+i +ai05 - 06,---- (2.42) 

The equations (2.34), (2.36) and (2.41) give rise to 



which indicates that the Hamilton- Jacobi equation is separable. Replacing i/^ by -q-^ 
and interpreting the Pi as integration constants, one gets the generating function S of 
the canonical transformation from (2.43) 



A^fe 



W{X) 



N+2 



The linearizing coordinates are then 

dS iV^ f^'^ X' 



The linear flow induced by (2.10) is then given by (using (2.42)) 

dFi 

Qi = 'li + x——=-ii-x5i^N+i, i = 0,1,...,A^+1, (2.46) 

where 7^ are arbitrary constants. Combining the equation (2.45) with the equation 
(2.46) leads to the Jacobi inversion problem for the FDIHS (2.10) 

-V/ =dX = -fi-x6^N+l, z = 0,l,...,A^+l. (2.47) 

Since V'li ? ^i 7 92 defined by (2.39) are the symmetric functions of Hk,k = 1,...,N + 2 
by using the standard Jacobi inversion technique [33], they can be solved in terms of 
Riemann theta functions from (2.47). After having ipij,qi,q2, the ip2j,PiiP2 can be 
found by using (2.10a). In this way the solution to (2.10) is obtained. 
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III. The first family of QBH systems. 

In the foUowing sections, by using the method described in the previous section, we 
wiU present QBH representation for some families of FDIHSs given in [24], and prove 
the equivalence of two sets of separated variables. 

A. The first family of FDIHSs. 

We first recall the constrained flows of the hierarchy of nonlinear evolution equations 
(NLEE) associated with the following polynomial second order spectral problem [31] 



ij, = U{u,X)^, U{u,X) 







1 



ET=o^^y 



i; 



(3.1) 



where Um = —l,u = (um-i, •••,tio)^- The adjoint representation (2.2) of (3.1) yields 



flo 



O-m — bo — ... — bm-1 — 0, bm — 1, ^m+l — ^"Wm-l, 



1 1 

flm+l — —^Um-l,xi Co — 1, Ci — —-Um-1, • • • , 



and in general 



/ bk+m \ 



\ bk+1 J 



I bk+m-l \ 



\ h J 



(3.2a) 



where 



1 1 

dk = —-^bk,x, Ck = —-bk,xx — 2_^'^i^k+i, k = l,2, 

1=0 



f Lm-1 Ltji-2 

1 

1 



V 



1 9 1 1 

Lo = -D^ + uo - -D-\o,^ 



1 



Ui D u 



Li Lo\ 




1 0/ 

1, 



(3.2b) 



i = 1, . . . , m — 1. 



The hierarchy of NLEEs associated with (3.1) can be written as an infinite- dimensional 
Hamiltonian system 



«t. 



( Um-l \ 



\ «0 /. 



/ bn+m \ 



J 



J- 



\ bn+1 J 



Su 



n = l,2,..., 



(3.3) 
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where the Hamiltonian H„, and the Hamiltonian operator J are defined by 



J 



^ 







2D \ 

2D Jm-l 

•J-m — l '-'m — 2 



\2D Jm-l 

Ji^ -Ui^^-2uiD, i = 0, l,...,m- 1, H, 
Under zero boundary condition we have 



Ji Jo / 

2 



m — 2n — 2 



y iuih 



'n+i+l- 



i=l 



6u 



(A-v?,A-v^...,^?r, 



SX _ dX 

du Su 



(3.4) 



Similarly, the constrained flows of the NLEEs (3.3) are deflned by [24] 



m — l 



*l,x = *2, *2,x = A"^*l - J2 ^^^'^^ 



i=0 



SH, 1 A SX 



E 



1 



5u 2 ^^-^ 5u 



h+ij 



0. 



(3.5a) 



(3.5b) 



<^i,^i> / 



For / = m, (3.5b) leads to 



Ur 



-^ = E(-i) 



i=i 



-i J + 1 
2^' 



^ < A'i^i,^i > ... < A'^^1,^1 >, 



ii+...+Zj=fc— j 

fc = 1, . . . ,m, 



(3.6) 



where /i > 0, ...,lj > 0. By substituting (3.6) into (3.5a), the flrst constrained flow of 
(3.3) can be written as a canonical FDHS 



^1, 



dFo 



^ 



dFo 



2,x 



or 



Px = 0o\7 Fq, 



(3.7) 



where 



p = (vl/f,vl/^)^ ^0 



InxN 

-Inxn 



FAMILIES OF QUASI-BI-HAMILTONIAN SYSTEMS AND SEPARABILITY 17 

1 "" _1 

Fo=-<*2,*2>+$^( — )^ + ' J2 <A'^*i,*i>...<A'^+^^l,*l>. 

j=0 li + ...+lj + i=m-j 

The entries of the Lax matrix for (3.7) are given by [24] 

m k ^ 

C(A) = A"^ + ^A"^-'=^(--)-'' Yl < A''^i,^i >...< A'^^1,^1 > 

A:=l j = l h + ...+lj=k-j 

We have 

A(A)2 + i?(A)C(A) = X-' + J2 X^A"' ^3.9) 

i=l 
^ m k ^ 

F« = _[A- + ^ Ar' ^(-2^' E < A^^vEr,, vi>, > . . . < A'. vi>,, vi>, >]^2. 

k=i j=i ii+...+ij=k-j 

1 ,2 , 1 V^ (^H^2fc - V'lfc^2i)^ ■ 1 .. 

where F*^*'', i = 1, ..., A^, are independent integrals of motion for (3.7) and Fq = X]i=o -^ • 
It can be shown that the system (3.7) is integrable in the LiouviU's sense. The systems 
with TO = 1, 2, ... give rise to a family of FDIHSs which include the well-known Garnier 
system as the first member (m=l). This family of FDIHSs was first given in [34]. 

In order to find the QBH structure for (3.7), we need to consider the following 
modified polynomial second order spectral problem [31] 

^. = [/(„, AM, f/(.,A)=(_^£«^^^._, _^^J. i.= (^^), (3,10) 

where Vm = —i,v = (vq, ...,Vm-i)'^ ■ The equations (2.2) and (2.3) yield 

1 

ao — . . . — am-2 — bo — ... — bm-S — 0, bm-2 — li bm-l — -Vm-1, 

1 

am-l=Vo, Co = 1, Ci = --Vm-1,- ■ ■, 
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and in general 



I 2ak+i \^ 



V —bk+m-l ' 



L 



( 2afc N^ 
-bk 



\ —bk+m-2 I 



Cfc+i = aA;,x - ^t'i^fc+i-i, /c = l,2, 



i=\ 



where 



/ -2t)o + D 
1 



L = 







. 

L2 . 








\ 




1 

Lm—2 Lm-1 / 



vu 



VO 



' ^m — 1 / + 



J 



( 2a„ ^^ 



I 



J 



5Hn 

Su 



n = l,2, 



(3.11a) 



(3.11b) 



Lo = -D+-D-^voD, L, = -v, + -D-^v,D, z = l,...,m-l. 
The hierarchy of NLEEs associated with (3.10) is 



(3.12) 



\ ~bn+m-2 . 

where the Hamiltonian H„, and the Hamiltonian operator J are given by 



/iD 



J 



Jo 



■h 



Jrn-1 -2D 
\ -2D 



\ 

• Jm-l —2D 



/ 



Ji = Vi,j: + 2viD, z = 0,1, . . .,m- 1, H, 



m — 2n — 2 



[an,x - y^JVjbn+i-l] 



1=1 



Also we have 



f = {2<p,<p„<plX<Pl...,X—'<Plf. 
ou 



(3.13) 



The constrained flows of (3.12) are defined by 



m — 1 



$l,x = Vo^i + A$2, $2,x = (A"" ^ - XI ^*^' ^)*1 " ^0*2, 



(3.14a) 



i=l 
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N 



Sv 2^ 6v 
i=i 



/ 2az \ 
-hi 



+ 



For I = m — 1, (3.14b) leads to 



/ 2<$i,$2> \ 

< $i,$i > 



(3.14b) 



Vo = --< $1,$2 >, 



(3.15a) 



'^m—k 



E(-i) 



i=i 



-i i + 1 

2^' 



Zi + ...+Zj=fc— J 

A; = 1 , . . . , m — 1 . 



(3.15b) 



By substituting (3.15) into (3.14a), the first constrained fiow of NLEE (3.12) can be 
written as a canonical FDHS 



or 



where 



$ 



l,x 



dFo dFp 

a$2' ^'"^ 9$i 

Px = 6'o V -^0 7 



^T /rTnT 



p = ($^$^)^ ^0 



Inxn 

-Inxn 



Fo = - < A$2, $2 > -^ < $1, $2 >^ 

m ^ 

+ ^(--)-'' J]] < A'i$i,$i > ...< A'^$i,$i > 



i=i 



ii + ...+Zj=m— J 



(3.16a) 



(3.16b) 



B. The QBH structure for the family of FDIHD (3.7). 

It is known [31] that the gauge transformation between the spectral problems (3.1) 
and (3.10) is given by 

Ui = Vi, z = l,...,m-l, uo = -vq^-Vq, (3.17) 

which, together with (3.6) and (3.15), gives rise to the map relating (3.7) to (3.16), i.e. 
P^ M{P): 



^1 = $1, ^2 = A$2 - - < $1, $2 > $1- 



(3.18) 
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In fact the map M transforms the first equation and the second equation with an 
additive term , — c\l/i(c = Fq), in (3.7) into the corresponding equations in (3.16). Since 
the 9i constructed in the following is vahd for an arbitrary c, so we can take c = 0. The 
Jacobi M' of the map M takes the form 

M'{P)=( 1 ^ ^ ^^"/ 1^ ^T ^^TJ'^T]■ (3.19) 

Then the second compatible Poisson tensor for the vector field (3.7) is 

By a straightforward calculation, we have the following proposition. 
Proposition 3. The system (3. 7) possesses the QBH representation 

P, = eoVFo = -01 V El (3.21a) 

P 

where 

p = S(A)U=o = 1 - ^ < A"'*i, ^1 >, (3.21b) 

TV 

El = [AHX)+B{X)C{X)]\,^o = -E^"'^^^^ 

1 1 

= - < A"^^2, *2 > +t[< A"^*1, *2 >^ - < A"^*l, *l >< A-^^2, ^2 >] 

m ^ 

+ ^(--)^+^ J2 < A'i^i,*i > ... < A^^+i-i*i,*i > . (3.21c) 

j=0 li + ...+lj + i=m-j 

C. The Nijenhuis coordinates. 

In the same way as for (2.30)-(2.33), the eigenvalues of the Nijenhuis tensor /ii, ..., fx^ 
are defined by the roots of the equation 

f{X)=\XI-Ai\=0, (3.22a) 

which gives 
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Then one defines 



dS 



N 



E^ 



99j 



2j- 



9f^j f~[ dfij 



j = l,...,N. 



(3.22b) 



On the other hand, the generahzed eUiptic coordinates (/i,i>) are defined by means 
of the Lax matrix in the foUowing way [24]. The coordinates /ii, ...,jiN are introduced 
by the zeros of B{X): 



1 

b{x)=i+-y: 



^h _ R{X) 



2^^A-A, K(A)' 
where K{X) is defined by (2.35) and 

AT N 

i?(A)=n(A-/i.)=E/^^^'^"'' 



(3.23a) 



(3.23b) 



fe=i 



fe=0 



N N 

j=i i=i 

and the canonicaUy conjugate coordinates t^i, ..., P^r are defined by 



N 



2 ^ /xfc - A, 



t^k = -A{fik) = 7;Y1 -'''''7 ^ k^l,...,N. 



(3.23c) 



We have the following proposition. 



Proposition 4. The Nijenhuis coordinates {p,v) defined by (3.22) are exactly the same 
as the generalized elliptic coordinates {p>,v) defined by (3.23). The QBH vector field 
(3.21) is Pfaffian in the Nijenhuis coordinates. 

Proof. Similarly, we have by induction 

/Ar(A; Ai, ..., Xn) = \XI — All 



X-Xi + ^1pfl ^^llV^12 

|V'12V'11 A - A2 + i^i2 



|^12^17V 



^^liV^ll li'lNi'12 ■■■ X-XN + li'^N 

1 ^ 

(A - Ai)/Ar-i(A; A2, ..., Xn) + -^?i 1[{X- A^) = B{X)K{X), (3.24) 

fc=2 
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which shows that Hk = fik- It follows from (3.23a) that 

(3.25) 



Thus we have 


^?.^2^^, j^l,...,N. 




N r. 


/2i?(A,) ^ i^2jR{Xj) 


~ 2^^l^k- Xj 



, (3.26) 



which implies that Uk = I'k, since nu = Hk- Finally, it is found from (3.23a) that 

p = S(A)|a=o = 1 - ^ < A-'^i, ^1 >= — . 

2 CtAT 

This completes the proof. 

IV. The second family of QBH systems. 

For / = m + 1, it is found from (3.5b) [24] that 

U^-k = (-2)''«m-l 



k-2 [^-1 



i=0 j = l ;^_|_..._|_;^=fc_i_2j 

(4.1a) 

m — l 

Lo«^-i =< A— 1^1, ^1 >-J2L^< A^-^^1, ^1 >, (4.1b) 

where 

Z i.J . 

Denote 

_ 1 

Q. ^m — 1? P TJ'^m— l,x- 

o 

By substituting (4.1a), (3.5a) and (4.1b) become a canonical FDHS 

Px^^oV^i, (4.2a) 

where 

\^-'(Af+l)x(Ar+l) U 
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rrri + 2- 
m [ 2~ 



+ 5Z^' J2 f^-'J Y. <A'^*i,*i>...<A'^*i,*i>. (4.2b) 

The entries of the Lax matrix Q for (4.2) are of the form [24] 

^W-2,- EMf. B(A)^A. ,. E,^. (^.Sa) 

m+1 JV ,2 

(7(A)=$:A-+^-'=c.-i5:-^, (4.3b) 

fc=0 j = l -^ 

where 

^ A;-2 [^-] 

i=0 j=l li + ...+lj=k-i-2j 

k = l,...,m+l, (4.3c) 

g^+2+fc = -- < A'^^2,^2 >, /c = 0,l,.... (4.3d) 



Similarly, the equality 



^2(A) + S(A)C(A) = A"^+2 - Fi + ^ -— -, (4.4) 



. , A - A, 

1=1 



1 1 1 """^^ 



k=0 



1 V^ (^li'02A; - ll^lk'4^2i)'^ 



+ 4^1. A^^A^ ' ^ = 1'-'^' 

determines A^ + 1 independent integrals of motion Fq, F^^\ ^ = 1, ..., N^ for the FDHS 
(4.2). The systems (4.2) for m = 1,2, ..., give the second family of FDIHSs. By taking 
m — 1 (4.2) gives rises to the multidimensional Henon-Heiles system. The system (4.2) 
was also studied by a recurrence relation in [35], however no explicit expressions like 
(4.2b) and (4.3) were given in that paper. 

In the exactly the same way as we did in the previous section, we can obtain another 
FDHS from (3.14) for / = m, find the map relating this FDHS to the FDHS (4.2) and 
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finally, by using this map, obtain the second compatible Poisson tensor for the vector 
field for (4.2) 



-AJ B 



A -i*i\ „ _ f OnxN 1*2 

By a straightforward calculation, we can show the following proposition. 
Proposition 5. The system (4-2) possesses the QBH representation 

P. = eoVFi = -e, V El (4.6) 

p 



where 



p = B{X)\x=o = \^-\< A"'*i, *i >, (4.7a) 

N 



El = [A\\) + S(A)C(A)]|a=o = -Fi-Y, K'F 



i=l 



2p < A-i^i, ^2 > +]q < A"'*2, *2 > +V - {-\qr^^ 



r rn, + 2- 



5Z^' S^ ^^'^ I] <A'^*i,*i>...<A^^*i,*i> 

1=0 j = l h + ...+lj=m+2-i-2j 



r T7i-|-1 — i -\ 

^ ^ m+1 L 2 J 

+ -[<A-l^l,^2>'-<A-l^l,^l><A-1^2,*2>]--5^g' 5Z P^,J 

i=0 j=0 

X ^ < A^i^i,^i > ... < A'^^1,^1 >< A^^+1-^^1,^1 > . (4.7b) 

ll + ...+lj + i=m+l — i — 2j 

In the same way, pi, ...,fxis[+i in the Nijenhuis coordinates are defined by the roots 
of the equation 

fiX)=\XI-Ai\=0, (4. 

which gives 
Then one defines 



".-^-1:^.11-^^. i-i.....^-i- (-b) 
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On the other hand, the generahzed parabohc coordinates {fi,, u) are defined by means 
of the Lax matrix in the foUowing way [24]. The coordinates /ii, ..., JiN+i are introduced 
by the zeros of B{X): 

^W-A.l,.l|:3A_.m (,9a) 

where K{X) is defined by (2.35) and -R(A) by 

N+l N+1 

R{X)= J] (A - ^fc) = 5^ /3,A^+l-^ 

fc=l fc=0 

N N+l 

/3o = l, A = -5^^,,..., f3N+i = {-l)''+'l[flj, 
and the canonically conjugate coordinates ui, ..., i'n+i are defined by 

% = -A(/Xfc) = -2p+^V^^i^, k=l,...,N+l. (4.9b) 

We have the following proposition. 

Proposition 4. The Nijenhuis coordinates {lJi,v) defined by (4-8) are exactly the same 
as the generalized parabolic coordinates [p-.v] defined by (4-9). The QBH vector field 
(4-6) is Pfaffian in the Nijenhuis coordinates. 

Proof. In a similar way, we can show by induction that 

/(A) = B{\)K{\). (4.10) 

It follows from (4.9a) that 



-,2 _ r, R{^3 



V'h=2t777VT, i = l,-,iV, 



N 



OiN ^ XjK'{Xj) aN 



3 

Then it is similar to find that 

p = B{X)\x=o = Iq-1< A-'*i, *i >= ^^^ 

2 2 aN 

This completes the proof. 
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V. Concluding remarks. 

In the exactly same way as we did in the previous two sections, we can construct the 
third family of QBH systems from the constrained flows (3.5) for / = TO + 2,m= 1,2,.... 
The QBH system (2.28) is just the second member {m = 2) in the third family of QBH 
systems, and 9i and p given by (2.27) and (2.29a) are the second compatible Poisson 
tensor and the integrating factor for the third family of QBH systems. 

In general, the constrained flow (3.5) for I = m + k can be transformed into a 
FDIHS by introducing the Jacobi-Ostrogradsky coordinates. Under the map relating 
this FDIHS to that obtained from the modified constrained flow (3.14) for / = m + k — l, 
the image of the Poisson tensor ^o for the latter gives rise to the second compatible 
Poisson tensor 9i for the former. In this way, for each k we can obtain a family of 
QBH systems with m = 1,2, .... The results obtained in the previous sections suggest 
the following conjecture: each family of QBH systems {I = m + k,m = 1,2, ...) shares 
the same 9i and p for the QBH structure 

^oV^i = -^iV^i, 
P 

and, in general, by mean, of the Lax matrix Q = ( ^f,'l\ ''^^l, ) and the expression 

\C[\) -A{\)J 



m+2k N 



A'{x) + B{x)c{x)= J2 F^y + Yl 

i=0 1=1 

we have 



A -A,; 



AT 



p = i?(A)U=o, E^ = [A^X) + B{X)C{X)]\x=o = Fo - J^F^A-^. 

i=l 

For A; = 1, 2, .., we find an infinite number of families of QBH systems. Furthermore we 
can show in a similar way that the Nijenhuis coordinates introduced by the Nijenhuis 
tensor are exactly the same as the separated variables defined by means of the Lax 
matrix for the QBH system in the family, and each QBH vector field is Pfaffian in the 
Nijenhuis coordinates. 
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